Probing solution conformations of l-DOPA: Integration of experiment and simulation via vibrational optical activity.

l-DOPA plays a critical role as a precursor to dopamine and is a standard treatment for Parkinson's disease. Recent research has highlighted the potential therapeutic advantages of deuterated l-DOPA analogs having a longer biological half-life. For their spectroscopic characterization, the in-detail characterization of l-DOPA itself is necessary. This article presents a thorough examination of the vibrational spectra of l-DOPA, with a particular emphasis on chirally sensitive VOA techniques. We successfully obtained high-quality Raman and ROA spectra of l-DOPA in its cationic form, under low pH conditions, and at a high concentration of 100 mg/ml. These spectra cover a broad spectral range, allowing for precise comparisons with theoretical simulations. We also obtained IR and VCD spectra, but they faced limitations due to the narrow accessible spectral region. Exploration of l-DOPA's conformational landscape revealed its intrinsic flexibility, with multiple coexisting conformations. To characterize these conformations, we employed two methods: one involved potential energy surface scans with implicit solvation, and the other utilized molecular dynamics simulations with explicit solvation. Comparing ROA spectra from different conformer groups and applying spectral decomposition proved crucial in determining the correct conformer ratios. The use of explicit solvation significantly improved the quality of the final simulated spectral profiles. The accurate determination of conformer ratios, rather than solely relying on the number of averaged spectra, played a crucial role in simulation accuracy. In conclusion, our study offers valuable insights into the structure and conformational behavior of l-DOPA and represents a valuable resource for subsequent spectroscopic studies of its deuterated analogs.

Structure of Zinc and Nickel Histidine Complexes in Solution Revealed by Molecular Dynamics and Raman Optical Activity.

The histidine residue has an exceptional affinity for metals, but solution structure of its complexes are difficult to study. For zinc and nickel complexes, Raman and Raman optical activity (ROA) spectroscopy methods to investigate the link between spectral shapes and the geometry were used. The spectra were recorded and interpreted on the basis of ionic equilibria, molecular dynamics, ab initio molecular dynamics, and density functional theory. For zwitterionic histidine the dominant tautomer was determined by the decomposition of experimental spectra into calculated subspectra. An octahedral structure was found to prevail for the ZnHis2 complex in solution, in contrast to a tetrahedral arrangement in the crystal phase. The solution geometry of NiHis2 is more similar to the octahedral structure found by X-ray. The Raman and ROA structural determinations of metal complexes are dependent on extensive computations, but reveal unique information about the studied systems.

Effects of sulfation and the environment on the structure of chondroitin sulfate studied via Raman optical activity.

Glycosaminoglycans are linear carbohydrate polymers with essential roles in many biological processes. Chondroitin sulfate (CS) is one of them, omnipresent in living organisms as an important structural component of cartilage. It provides much of its resistance to compression. Despite its biological importance, little is still known about the relation of the CS structure to chemical composition and interaction with the environment. We therefore measured Raman and Raman optical activity (ROA) spectra of five CS samples of different biological origin and variously sulfated CS building blocks (GlcA, GalNAc, and basic disaccharide units) in a wide frequency range between 200 cm-1 and 1800 cm-1 and analyzed them with respect to specific structure marker bands. We show that ROA spectroscopy is sensitive to the conformational stability and rigidity of pyranose rings of saccharides, the orientation of sugar hydroxyl groups and the secondary structure of the CS's backbone. The CS secondary structure has been found to be quite stable, with a minor variation as a reaction to physicochemical parameters (concentration, pH, temperature, and the presence of cations). Larger changes were observed under chemical changes (sulfation) of the CS chain. ROA spectroscopy thus exhibited useful potential to study the structure of similar biopolymers.

Absolute Configuration Determination of a Taxol Precursor Based on Raman Optical Activity Spectra.

Enantiomeric purity of drugs is essential for their biological activity. In the present study, we investigate the performance of Raman optical activity (ROA) spectroscopy in distinguishing four possible stereoisomers of the synthetic precursor used for the production of taxol from baccatin III. Taxol is one of the best-selling medicaments used in the treatment of ovarian, lung, and breast cancers and Kaposi's sarcoma. In a low yield, it may be isolated from the bark of the Pacific yew tree (Taxus brevifolia); however, its industrial production is largely dependent on the precursor. It is shown that for reliable comparison of the experimental and computed Raman and ROA intensities a large number of conformers had to be averaged, to properly account for molecular flexibility in solution. In addition, if combined with the density functional theory computations, ROA spectra provide convenient and economic means of absolute configuration determination.

Organization of the MADS box from human SRF revealed by tyrosine perturbation.

MADS box family transcription factors are involved in signal transduction and development control through DNA specific sequence recognition. The DNA binding domain of these proteins contains a conservative 55-60 amino acid sequence which defines the membership of this large family. Here we present a thorough study of the MADS segment of serum response factor (MADS(SRF)). Fluorescence, UV-absorption, and Raman spectroscopy studies were performed in order to disclose its behavior and basic functional properties in an aqueous environment. The secondary structure of MADS(SRF) estimated by analysis of Raman spectra and supported by CD has revealed only the C-terminal part as homologous with those of free core-SRF, while the N-terminal part has lost the stable α-helical structure found in both the free core-SRF and its specific complex with DNA. The three tyrosine residues of the MADS(SRF) were used as spectroscopic inner probes. The effect of environmental conditions, especially pH variations and addition of variously charged quenchers, on their spectra was examined. Two-component fluorescence quenching was revealed using factor analysis and corresponding Stern-Volmer constants determined. Factor analysis of absorbance and fluorescence pH titration led to determination of three dissociation constants pKa1 = 6.4 ± 0.2, pKa2 = 7.3 ± 0.2, and pKa3 = 9.6 ± 0.6. Critical comparison of all experiments identified the deprotonation of His193 hydrogen bonded to Tyr195 as a candidate for pKa1 (and that of Tyr158 as a candidate for pKa2). Within MADS(SRF), His193 is a key intermediary between the N-terminal primary DNA binding element and the hydrophobic C-terminal protein dimerization element.

Electronic circular dichroism of the chiral rigid tricyclic dilactam with nonplanar tertiary amide groups.

Electronic circular dichroism (ECD) of the spirocyclic dilactam 5,8-diazatricyclo[6,3,0,0(1,5)]undecane-4,9-dione has been measured in the extended wavelength range (170-260 nm) utilizing far-UV CD instrumentation including synchrotron radiation light source. The data of this model of two nonplanar tertiary amide groups interacting within the rigid chiral environment provided new information particularly about the shorter wavelength π-π* transition region below 190 nm. The interpretation using TDDFT calculations confirmed that effects of amide nonplanarity follow our previous observations on monolactams as far as amide n-π* transitions are concerned. ECD band in the n-π* transition region of the nonplanar diamide exhibits an identical bathochromic shift and its sign remains tied to the sense of nonplanar deformation in the same way. As far as n-π* transitions are concerned amide nonplanarity acts as a local phenomenon independently reflecting sum properties of single amide groups. On the other hand, CD bands associated with π-π* transitions (found between ∼170 to 210 nm) form an exciton-like couplet with the sign pattern determined by mutual orientation of the associated electric transition moments. This sign pattern follows predictions pertaining to a coupled oscillator. The influence of amide nonplanarity on π-π* transitions is only minor and concentrates into the shorter wavelength lobe of the π-π* couplet. The detailed analysis of experimental ECD with the aid of TDDFT calculations shows that there is only little interaction between effects of inherent chirality caused by nonplanarity of amide groups and amide-amide coupling. Consequently these two effects can be studied nearly independently using ECD. In addition, the calculations indicate that participation of other type of transitions (n-σ*, π-σ* or Rydberg type transitions) is only minor and is concentrated below 180 nm.

Protonation effect of tyrosine in a segment of the SRF transcription factor: a combined optical spectroscopy, molecular dynamics, and density functional theory calculation study.

The high sensitivity to pH of a short segment (an octamer) of serum response factor (SRF), an important member of the MADS box family of transcription factors, was investigated by Raman scattering, infrared and circular dichroism spectroscopies. Molecular dynamics (MD) and density functional theory (DFT) calculations enabled interpretation of spectral changes in close detail. Although there was a negligible difference between spectra in acidic and neutral environments, the spectrum in basic pH was substantially different. The major changes were attributed to the deprotonation of tyrosine. The secondary structure of the SRF octamer fragment was estimated experimentally as well as predicted theoretically by MD. All techniques proved that it exists in a dynamical equilibrium among several conformations mostly close to ß turn, unordered conformations, and extended structure, in contrast to the stable secondary structure it possesses as a part of SRF. Generally, this approach represents a useful tool for the study of various short oligopeptides.

Nonplanar tertiary amides in rigid chiral tricyclic dilactams. Peptide group distortions and vibrational optical activity.

We investigate amide nonplanarity in vibrational optical activity (VOA) spectra of tricyclic spirodilactams 5,8-diazatricyclo[6,3,0,0(1,5)]undecan-4,9-dione (I) and its 6,6',7,7'-tetradeuterio derivative (II). These rigid molecules constrain amide groups to nonplanar geometries with twisted pyramidal arrangements of bonds to amide nitrogen atoms. We have collected a full range vibrational circular dichroism (VCD) and Raman optical activity (ROA) spectra including signals of C-H and C-D stretching vibrations. We report normal-mode analysis and a comparison of calculated to experimental VCD and ROA. The data provide band-to-band assignment and offer a possibility to evaluate roles of constrained nonplanar tertiary amide groups and rigid chiral skeletons. Nonplanarity shows as single-signed VCD and ROA amide I signals, prevailing the couplets expected to arise from the amide-amide interaction. Amide-amide coupling dominates amide II (mainly C'-N stretching, modified in tertiary amides by the absence of a N-H bond) transitions (strong couplet in VCD, no significant ROA) probably due to the close proximity of amide nitrogen atoms. At lower wavenumbers, ROA spectra exhibit another likely manifestation of amide nonplanarity, showing signals of amide V (δ(oop)(N-C) at ~570 cm(-1)) and amide VI (δ(oop)(C'═O) at ~700 cm(-1) and ~650 cm(-1)) vibrations.

CH Stretching Region: Computational Modeling of Vibrational Optical Activity.

Most organic compounds provide vibrational spectra within the CH stretching region, yet the signal is difficult to interpret because of multiple difficulties in experiment and modeling. To better understand various factors involved, the ability of several harmonic and anharmonic computational approaches to describe these vibrations was explored for α-pinene, fenchone, and camphor as test compounds. Raman, Raman optical activity (ROA), infrared absorption (IR), and vibrational circular dichroism (VCD) spectra were measured and compared to quantum chemical computations. Surprisingly, the harmonic vibrational approach reasonably well reproduced the measured spectral patterns, including the vibrational optical activity (VOA). The CH stretching, however, appeared to be more sensitive to the basis set and solvent variations than lower-frequency vibrations. For a higher accuracy in frequencies and spectral shapes, anharmonic corrections were necessary. Accurate harmonic and anharmonic force fields were obtained with the mPW2PLYP double-hybrid functional. A limited vibrational configuration interaction (LVCI) where the CH stretching motion was decoupled from other vibrations provided the best simulated spectra. A balanced harmonic oscillator basis set had to be used, containing also states indirectly interacting with fundamental vibrations. A simpler second-order perturbational approach (PT2) appeared less useful. The modeling provided unprecedented agreement with experimental vibrational frequencies; spectral shapes were reproduced less faithfully. The possibility of ab initio interpretation of the CH spectral region for relatively large molecules further broadens the application span of vibrational spectroscopy.

Tracking of the polyproline folding by density functional computations and Raman optical activity spectra.

Polyprolines offer many opportunities to study factors influencing peptide and protein folding and structure. Longer chains can adopt two well-defined forms (PPI and PPII), but shorter peptides are quite flexible. To understand in detail the dependence of the secondary structure on the length and the interplay between the side chain and main chain conformation, zwitterionic (Pro)(N) models (with N = 2, 3, 4, 6, 9, 12 and longer inhomogeneous chains) were studied by a combination of the Raman and Raman optical activity (ROA) spectroscopy with the density functional theory (DFT). Potential surfaces were systematically explored for the shorter oligoprolines, and Boltzmann conformational ratios were obtained both for the main chain and the proline ring puckering. The predictions were verified by comparison of the experimental and simulated ROA spectra. The conformer ratios extracted from a decomposition of the experimental ROA into scaled computed spectra well reproduced Boltzmann populations calculated from relative energies. For example, an "A" puckering of the proline ring was found prevalent, relatively independent of the length, whereas the cis-amide backbone form adopted by shorter peptides rapidly disappeared for N > 4. The results are consistent with previous NMR and vibrational circular dichroism (VCD) data. Delocalized exciton vibrations along the peptide chain often enhance the ROA signal, and can thus be used to indicate a longer regular peptide structure. The ROA technique appeared to be very sensitive to the ring puckering; less distinct spectral features were produced by changes in the main chain geometry.

Photoinduced processes in hydrogen bonded system: photodissociation of imidazole clusters.

The photodissociation of imidazole in hydrogen bonded clusters has been studied at photodissociation wavelengths 243 and 193 nm. Imidazole clusters of different mean cluster sizes n approximately 3 and 6 have been produced in expansions with He and Ar carrier gases, and the mean cluster sizes were determined by mass spectrometric and crossed beam scattering experiments. Simultaneously, the (C(3)N(2)H(4))(n) clusters were studied by ab initio calculations for n up to 4 molecules, confirming the hydrogen bond N-H...N motif in the clusters. The measured H-fragment kinetic energy distribution spectra exhibit a bimodal character similar to the KEDs found for the bare molecule. (1) At 243 nm the fast H-atoms originate from the direct dissociation process on the repulsive pi sigma* state, and the slow component results from the dynamics populating the vibrationally hot ground state via an S(1)/S(0) conical intersection. In the clusters the contribution of the slow component increases with the cluster size. The slow component is also dominant at the shorter wavelength of 193 nm, where the dynamics starts with the excitation of pi pi* state. It is shown that the slow component in our experiment is a product of subsequent two-photon absorption. We have proposed different mechanisms how the observed enhanced internal conversion can be rationalized. The increased stability with respect to the H-fragment dissociation in clusters can be caused either by hydrogen transfer in the N-H...N bond or by closing the pi sigma* dissociation channel as in the case of pyrrole clusters.

Fragmentation dynamics of size-selected pyrrole clusters prepared by electron impact ionization: forming a solvated dimer ion core.

Structure and dynamics of size-selected charged pyrrole clusters have been studied by means of molecular beam scattering experiments and ab initio calculations. Small neutral Pyn clusters were produced in Py/He mixture expansions, and the scattering experiment with a secondary beam of He-atoms was exploited to select the neutral clusters of different sizes. The complete size-selected fragmentation patterns for the neutral dimer to the tetramer after an electron impact ionization at 70 eV were obtained from the measurements of the angular and velocity distributions at different fragment masses. All the investigated cluster sizes decay mainly to the monomer ions Py+1 (from 60 to 80% of the corresponding neutral size) and to the dimer ion Py+2 (20-30%). The trimer ions Py+3 are generated to less than 10% from the neutral trimer and tetramer. To explain the observed results, we have calculated the structures and energetics of pyrrole clusters up to the trimer for the neutral and the ionic state using DFT and PMP2 methods. The ab initio calculations show that ionized pyrrole clusters are formed with a dimeric core that is solvated by neutral pyrrole molecules. In addition, the ground and ionic state of Py-Ar complexes were calculated at CCSD(T) level with extended basis in relevance to the mixed clusters produced in supersonic expansions of Py seeded in Ar. The calculated dissociation energies of the Py-Ar and (Py-Ar)+ complexes indicate that Ar atoms are able to rapidly evaporate after ionization. The combined analysis of the fragmentation probabilities, and calculations allowed us to estimate the distribution of energy deposited in the clusters after the electron impact, which peaks above 1 eV and has a tail up to 5 eV.

Experimental and theoretical study of the pyrrole cluster photochemistry: closing the pisigma* dissociation pathway by complexation.

Photolysis of size selected pyrrole clusters has been investigated and compared to the photolysis of an isolated pyrrole molecule. Experimentally, size distributions of different mean cluster sizes (n=3 and n>>5) have been prepared in supersonic expansions and the clusters were photolyzed at 243 and 193 nm. The kinetic energy distributions of the H photofragments have been measured. The distributions exhibit a bimodal character with fast and slow H-fragment peaks similar to the spectra of the bare molecule. However, with increasing cluster size the slow component gains intensity with respect to the fast one. A similar effect is observed with increasing the excitation energy from 243 to 193 nm. Theoretical calculations at the CASSCF/CASPT2 level have been performed for bare and complexed pyrroles (pyrrole is complexed with an argon atom and with another pyrrole unit). Combination of theoretical and experimental approaches leads to the conclusion that the direct dissociative pathway along the pisigma* potential energy surface in the N-H stretch coordinate is closed by the presence of the solvent molecule. This pathway is an important channel leading to the fast H atoms in the dissociation of the bare molecule. The solvent molecule influences significantly the electronic structure in the Rydberg-type pisigma* state while it has little influence on the valence states. The slow channel is mostly populated by the out-of-plane deformation mode which is also not influenced by solvation. We have also studied other possible reaction channels in pyrrole clusters (hydrogen transfer, dimerization). The present study shows that more insight into the bulk behavior of biologically relevant molecules can be gained from cluster studies.